Geometric Phase in Entangled Systems: A Single-Neutron Interferometer Experiment 
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The influence of the geometric phase on a Bell measurement, as proposed by Bertlmann et al. 
in [Phys. Rev. A 69, 032112 (2004)], and expressed by the Clauser-Horne-Shimony-Holt (CHSH) 
inequality, has been observed for a spin-path entangled neutron state in an interferometric setup. 
It is experimentally demonstrated that the effect of geometric phase can be balanced by a change 
in Bell angles. The geometric phase is acquired during a time dependent interaction with two 
radio-frequency (rf) fields. Two schemes, polar and azimuthal adjustment of the Bell angles, are 
realized and analyzed in detail. The former scheme, yields a sinusoidal oscillation of the correlation 
function S, dependent on the geometric phase, such that it varies in the range between 2 and 2y/2 
and, therefore, always exceeds the boundary value 2 between quantum mechanic and noncontextual 
theories. The latter scheme results in a constant, maximal violation of the Bell-like-CHSH inequality, 
where S remains 2\/2 for all settings of the geometric phase. 

PACS numbers: 03.75.Dg, 03.65. Vf, 03.65.Ud, 07.60.Ly, 42.50.Dv 



I. INTRODUCTION 

Since the famous 1935 Einstein-Podolsky- Rosen (EPR) 
gedanken experiment [l[ much attention has been de- 
voted to quantum entanglement [2j , which is among the 
most striking peculiarities in quantum mechanics. In 
1964, J. S. Bell Q introduced inequalities for certain cor- 
relations, which hold for the predictions of any hidden- 
variable theory applied Q. However, a dedicated ex- 
periment was not feasible at the time. Five years later 
Clauser, Home, Shimony and Holt (CHSH) reformulated 
Bell's inequality (BI) pertinent for the first practical 
test of the EPR claim [5J] . Thereafter polarization mea- 
surements with correlated photon pairs produced by 
atomic cascade 0, HJ and parametric down-conversion of 
lasers [j| [13, [HI , demonstrated a violation of the CHSH 
inequality. Up to date several systems [l^ [H, 0, [Hj] 
have been examined, including neutrons [161 ]. 

EPR experiments are designed in order to test local 
hidden variable theories (LHVTs) thereby exploiting the 
concept of nonlocality. LHVTs are a subset of a larger 
class of hidden-variable theories known as the noncon- 
textual hidden- variable theories (NCHVTs). Noncontex- 
tuality implies that the value of a measurement is in- 
dependent of the experimental context, i.e. of previous 
or simultaneous measurements [13, EH ■ Noncontextual- 
ity is a more stringent demand than locality because it 
requires mutual independence of the results for commut- 
ing observables even if there is no spacelike separation 
[19j |. First tests of quantum contextuality, based on the 
Kochen— Specker theorem [20| |. have been recently pro- 
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pose d 12 ll [22| and performed successfully using trapped 
ions |23| and neutrons 0, Hfl . 

In the case of neutrons, entanglement is not achieved 
between different particles, but between different degrees 
of freedom. Since the observables of one Hilbert spaces 
(HS) commute with observables of a different HS, the sin- 
gle neutron system is suitable for studying NCHVTs. Us- 
ing neutron interferometry p6l . [27j | single-particle entan- 
glement, between the spinor and the spatial part of the 
neutron wave function [161 ]. as well as full tomographic 
state analysis [28] , have already been accomplished. In 
addition creation of a triply entangled single neutron 
state [2t|, by applying a coherent manipulation method 
of a neutron's energy degree of freedom, has been demon- 
strated. 

The total phase acquired during an evolution of a quan- 
tal system generally consists of two components: the 
usual dynamical phase — 1/h J H(t)dt, which depends on 
the dynamical properties, like energy or time, and a ge- 
ometric phase 7, which is, considering a spin | system, 
minus half the solid angle (f2/2) of the curve traced out in 
ray space. The peculiarity of this phase, first discovered 
by M. V. Berry in 1984 (30|, lies in the fact that it does 
not depend on the dynamics of the system, but purely 
on the evolution path of the state in parameter space. 
From its first verification, for photons in 1986 [3l| and 
later for neutrons [32|. generalizations such as non adi- 
abatic [33j . noncyclic [34| . including the Pancharatnam 
relative phase [33, off-diagonal evolutions l36l. l37l . l38j |. 
as well as the mixed state case [H, 0, |4l], |42|. l43j]. have 
been established. 

The geometric phase in a single-particle system has 
been studied widely over the past two and a half decades. 
Nevertheless its effect on entangled quantum systems is 
less investigated. The geometric phase is an excellent 
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candidate to be utilized for logic gate operations in quan- 
tum communication [44[, due to its robustness against 
noise, which has been tested recently using superconduct- 
ing qubits [45| , and trapped polarized ultracold neutrons 
46] . Entanglement is the basis for quantum communica- 
tion and quantum information processing, therefore stud- 
ies on systems combing both quantum phenomena, the 
geometric phase and quantum entanglement, are of great 
importance [13, IH, Eg] . 

This article reports on an experimental confirmation 
for the violation of Bells inequality, relying on correla- 
tions between the spin and path degrees of freedom of 
a single neutron system, under the influence of the geo- 
metric phase. The geometric phase is generated in one of 
the complementary Hilbert spaces, in our case the spin 
subspace. We demonstrate in detail how the geometric 
phase affects the Bell angle settings, required for a vi- 
olation of a Bell inequality in the CHSH formalism, in 
a polarized neutron interferometric experiment. In sec- 
tion |TT] the theoretical framework, as developed in (HI, is 
briefly described for a spin-path entangled neutron state. 
Expectation values and Bell-like inequalities are defined 
and the concept of polar and azimuthal angle adjustment 
is introduced. Section IIIII explains the actual measure- 
ment process. It focuses on experimental issues such as 
state preparation, manipulation of geometric phase, joint 
measurements, as well as the experimental strategy. In 
the principal part data analysis and experimental results 
are presented. This is followed by sections [IV] and [V] con- 
sisting of discussion, conclusion and acknowledgments. 



II. THEORY 

A. Expectation Values 

First we want to clarify the notation, since numerous 
angles are due to appear in this article. Angles denoted as 
a are associated with path, and angles denoted as (3 with 
spin subspace. The ' symbol is used to distinguish differ- 
ent measurement directions of one subspace, required for 
a CHSH-Bell measurement @ (e.g. a and a' represent 
the measurement directions for the path subspace). In- 
dex 1 denotes polar angles, whereas index 2 is identified 
with azimuthal angles (e.g. Pi and (3[ are polar angles 
of the spin subspace). Finally, the - 1 symbol is used for 
adding it to an angle (e.g. — ct\ + tt). 

Following the notation given in [48] , in our experiment 
the neutron's wavefunction is defined in a tensor product 
of two Hilbert spaces: One Hilbert space is spanned by 
two possible paths in the interferometer given by |I), |II), 
and the other by spin-up and spin-down eigenstates, de- 
noted as l - ^) and referred to a quantization axis in- 
duced by a static magnetic field. Interacting with a time 
dependent magnetic field, the entangled Bell state ac- 
quires a geometric phase 7 tied to the evolution within 



the spin subspace [48| 

|* Bell ( 7 )) = -L(]I)®|^) + |II)®e^)). (I) 

As in common Bell experiments a joint measurement for 
spin and path is performed, thereby applying the projec- 
tion operators for the path 



Pl( a ) = \±a)(±a\, (2) 



with 



\ +a ) = cos^|I) +e M2 sin^|II) 

|_ a) = - sin ^L|I> cos^|II), (3) 

where ot\ denotes the polar angle and oti the azimuthal 
angle, and, for the spin subspace, 



with 



PUP) = |±/3)(±/3|, 



|+/3) = cos^|ft) + e^sm^| 



(4) 



-0) 



. 01 



Mn-^)+e^cos^>. (5) 



Introducing the observables 

A? (a) = Pl(a)-Pl{a) 

B s ((3) = F»(p) - P s _(/3) (6) 

one can define an expectation value for a joint measure- 
ment of spin and path along the directions a and j3 

E(a,)3) = (tf|iP(a)<8> B S ((3)\V) 
= — cos cti cos /3i — cos(a 2 — 02 + 7) sin a\ sin (3% (7) 
= - cos(ai - Pi) for (a 2 - P2) = —7. 



B. Bell-like Inequalities 

Next, a Bell-like inequality in CHSH-formalism [5( is 
introduced, consisting of four expectation values with the 
associated directions a, a' and (3, (3 1 for joint measure- 
ments of spin and path, respectively 

S(a,a',t3,P',y) = \E(a,(3) - E{a, (3') + E(a' , (3) 
+ E(ac',f3')\ 

= — sinai (cos(a2 — 02 — 7) sin/?i 

— cos(tt2 — 02l) s ^ n 0[) — cosai(cos/3i — cos/3i) 

— sina' 1 (cos(a 2 — 02 — 7) sm@i + cos(a' 2 — 0' 2 — 7) sm(3[) 

— cos a[ (cos Pi + cos P[ 



(8) 
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FIG. 1: The experimental apparatus for joint measurement of spinor and path degrees of freedom with respect to the geometric 
phase. The incident neutron beam is polarized by a magnetic field prism. The spin state acquires a geometric phase 7 during 
the interaction with the two rf-fields and is flipped twice. The beam block is required for measurements solely in one path (± z 
direction of the path measurement). Finally, the spin is rotated by an angle S (in the x, z plane), by a dc-spin turner, for a 
polarization analysis and count rate detection. The Bloch-sphere description includes the measurement settings of a. and /3(S), 
determining the projection operators, used for joint measurement of spin and path, a is tuned by a combination of the phase 
shifter (%) and the beam block, and (3 is adjusted by the angle 8. 



The boundary of Eq. (|8|) is given by the value 2 for any 
NCHVT [H(J. Without loss of generality one angle can 
be eliminated by setting, e.g., a. = {a\ = (±2 = 0), 
which gives 



S(a',/3,/3', 7 ) 



x'l ^cos(o/ 2 — fa — 7) sin/3i 



COStti 



+ cos(a' 2 - 02 ~ 7) sin pfj 
x (cos 0i + cos j3[ ) — cos j3i + cos (3[ 



(9) 



Keeping the polar angles a! x , Pi and /3[ constant at the 
usual Bell angles a{ = f , (3 X = ^ , fl[ = ^ (and az- 

P' 2 = 0) reduces S 



imuthal parts fixed at a' 2 — P2 
to 



5( 7 ) = \-V2-V2 cos 



7 



(10) 



where the familiar maximum value of 2y/2 is reached for 
7 = 0. For 7 = 7r the value of S approaches zero. 



1. Polar Angle Adjustment 



The polar Bell angles (3\, P[ and a[ («i — 0), yielding 
a maximum S*-value, can be determined, with respect to 
the geometric phase 7, by calculating the partial deriva- 
tives (the extremum condition) of S in Eq. lfTTT) : 

dS 

—— = sin Pi + cos a[ sin Pi — cos 7 sin a[ cos Pi = 
dpi 

dS 

—— = ~ sin P[ + cos a[ sin P[ — cos 7 sin a\ cos P[ = 
opi 

dS 



- = sin a[ (cos Pi + cos P[ ) 



(12) 



— cos 7 cos a[ (sin Pi + sin P[ ) = . 

The solutions are given by 

Pi = arctan(cos7) (13a) 
P'l = tt-/?i (13b) 

n 1 = — . (13c) 
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With these angles the maximal S decreases for 7:0^ 
\ and touches at 7 = ^ even the limit of the CHSH 
inequality S — 2. 



Here we consider the case when the azimuthal angles 
are kept constant, e.g., a' 2 — P2 = P'2 = i a 2 = 0), 
denoted as 



S{a'i,Pi,P'i n ) 



sin a\ ^cos 7 sin Pi + cos 7 sin P[ 
cos ai (cos Pi + cos P[ ) (11) 
cos /9i + cos . 



2. Azimuthal Angle Adjustment 

Next we discuss the situation where the standard max- 
imal value S = 2a/2 can be achieved by keeping the polar 
angles a[ , Pi and P[ constant at the Bell angles a[ = ? , 
Pi = j, Pi = (ai = 0), while the azimuthal parts, 
a' 2 , P2 and P' 2 {u2 = 0), are varied. The corresponding S 
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FIG. 2: Schematic representation of the spinor evolution, 
within the two successive rf flippers, on the Bloch sphere. 
The geometric phase 7 is given by minus half of the solid an- 
gle fi, traced out by the state vector, depending on the phases 
<f>i and 0n of the oscillating magnetic fields in the rf flippers. 



function is denoted as 



V2- 



— [cos(a 2 



ft -7) 



+ C0S(CV2 -0' 2 - 7)) 



(14) 



The maximum value 2y/2 is reached for 
02 = 0' 2 , and 

"2-^2 = 7( mod7r )- 



(15a) 
(15b) 



For convenience 02 = is chosen. 

The conditions expressed in Eq. (|13ap . Eq. (|13b|) and 
Eq. (|15bj) (see Bloch spheres in Fig.[l} are experimentally 
realized using the spin turner device and the neutron 
interferometer depicted in Fig.[TJ 



ter and a subsequent spin flip process in one sub beam: 
Behind the beam splitter (first plate of the IFM) the neu- 
tron's wave function is found in a coherent superposition 
of 1 1) and I II), and only the spin in |II) is flipped by the 
first rf- flipper within the interferometer (see Fig.[T]). 

The entangled state which emerges from a coherent 
superposition of |I) and |II) is expressed as I^Bcii) = 

75 (l 1 ) ® lit) + l n ) ® ei0I l4)), where the state vector of 
the neutron acquires a phase <pi during the interaction 
with the oscillating field, given by BW = B^ ] cos(ujt + 
(pi) ■ y (for a more detailed description of the generation 
of I^Beii) see Hi). 

B. Manipulation of Geometric and Dynamical 
Phases 

Within the rf-flipper placed inside the interferometer 
(path II), the neutron spin traces out a semi-great circle 
from |ff) to |4) on the Bloch sphere and returns to its 
initial state |ff) when passing the second rf-flipper (see 
Fig. The two semi-great circles enclose an angle <f>\ — 
4>\i (orange slice), and hence a solid angle £1 = 2(<f>j — (f>u). 
The solid angle Q yields a pure geometric phase 7 = 
— fi/2 as in (Hll l52|. Since we set <ftu = the geometric 
phase is given by 7 = (fn and the state is represented by 



IWt)> = ^(|i> 



in) 



(16) 



In our experiment the |ff) eigenstate (in path I and II) 
acquires dynamical phase as it precesses about the mag- 
netic guide field in + z-direction. After a spin-flip (only in 
path II) the |JJ.) eigenstate still gains another dynamical 
phase but of opposite sign compared to the situation be- 
fore the spin-flip. These phases manifest in a dynamical 
phase offset, which remains constant during the complete 
measurement procedure. 



III. DESCRIPTION OF THE EXPERIMENT 
A. State Preparation 

The preparation of entanglement between spatial and 
the spinor degrees of freedom is achieved by a beam split- 



C. Joint Measurements 

Experimentally, the probabilities of joint (projective) 
measurements are proportional to the following count 
rates 



N ++ (a,8) =JV++(a,C9i,0)) cx (* Exp ( 7 )|pP(a) ® 0)|* Exp ( 7 )) 

N + -(a,P)=N ++ (a,(p 1 +Tr,0))=N ++ (a,(rt,0)) « (*Exp(7)|£?(a) ® HiPi* °)l*E*p(7)> 

N- + (a,/3) =iV ++ ((a 1 +7r,a 2 ),(/3 1 ,0)) = N ++ ((a± ,a 2 ), (0i,O)) cx (tf Exp ( 7 )|pP(cvr, a 2 ) ® 0)|^ Exp (7)) 

N^(a,f3)=N ++ ((ai,a 2 ),(pt,0)) cx (* Exp ( 7 )|pP(c^, a 2 ) ® P s +(0t , 0)|* Exp ( 7 )). 



The expectation value of a joint measurement of A v (a) is experimentally determined from the count rates 
and B s (d) 

E(a,/3) = 

E( a ,P) = (*( 7 )\AP(a)®B a (3)\*( 1 )) (18) N ++ (a,l3)-N + -(a,l3)-N- + (ct,P)+N—(a,P) (19) 

N ++ (a,/3)+N + _(a,/3) + N- + (a,/3)+N—(a,l3) ' 
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With these expectation values S is defined by 

S = E(a,0) - E(a,/3') + E{a',0) + E(a',/3'). (20) 



D. Experimental setup 

The experiment was carried out at the neutron 
interferometer instrument S18 at the high-flux reac- 
tor of the Institute Laue-Langevin (ILL) in Grenoble, 
France. A sketch of the setup is depicted in Fig.Q] 
A monochromatic beam, with mean wavelength A = 
1.91 A(AA/A ~ 0.02) and 5x5 mm beam cross-section, 
is polarized by a bi-refringent magnetic field prism in 
z-direction [531 ]. Due to the angular separation at the 
deflection, the interferometer is adjusted so that only the 
spin-up component fulfills the Bragg condition at the first 
interferometer plate (beam splitter). 

As in our previous experiment (29| . the spin in path 
| II) is flipped by a rf-flipper, which requires two mag- 
netic fields: A static field Bq-% with Bq = Hbj/(2\fi\) and a 
perpendicular oscillating field = cos(wt+<()i) -y 

with amplitude B^' = 7rfi/(2r|/i|), where /i is the mag- 
netic moment of the neutron and r denotes the time 
the neutron is exposed to the rf-field. The oscillating 
field is produced by a water-cooled rf-coil with a length 
of 2 cm, operating at a frequency of lu/2tt = 58 kHz. 
The static field is provided by a uniform magnetic guide 
field Bq ~ 2mT, produced by a pair of water-cooled 
Helmholtz coils. 

The two sub-beams are recombined at the third crystal 
plate where |I) and |II) only differ by an adjustable phase 
factor e lx (path phase x is given by x = Nb c \D, with 
the thickness of the phase shifter plate D, the neutron 
wavelength A, the coherent scattering length b c and the 
particle density N in the phase shifter plate) . By rotating 
the plate, x can be varied systematically. This yields 
the well known intensity oscillations of the two beams 
emerging behind the interferometer. 

The O-beam passes the second rf-flipper, operating 
at ui/2tt — 29 kHz, which is half the frequency of 
the first rf-flipper. The oscillating field is denoted as 
cos ((u>/2)i + 4>ijj ■ y, and the strength of the 
guide field was tuned to about 1 mT in order to satisfy 
the frequency resonance condition. This flipper compen- 
sates the energy difference between the two spin compo- 
nents, by absorption and emission of photons of energy 
E = fuv/2 (see (H). 

Finally, the spin is rotated by an angle 5 (in the x, z 
plane) with a static field spin-turner, and analyzed due to 
the spin dependent reflection within a Co-Ti multi-layer 
supermirror along the z-direction. With this arrange- 
ment consisting of a dc-spin turner and a supermirror 
the spin can be analyzed along arbitrary directions in 
the x, z plane, determined by 5, which is measured from 
the z axis (see Fig.Q] and later front panel of Fig. 2] for 
intensity modulations due to ^-scans). 



E. Experimental Strategy 

1. Polar Angle Adjustment 

Projective measurements are performed on parallel 
planes defined by a 2 = a' 2 — 02 — (3' 2 — (see Fig.QJ. 
For the path measurement the directions are given by 
a : oi\ = 0,^2 = (Fig. [3]), and a.' : a[ — n/2, a' 2 = 
(Fig.®. 

The angle a, which corresponds to +z (and — z for 
cti = (Xi + it = Tc,a.2 = 0) is achieved by the use of 
a beam block which is inserted to stop beam II (I) in 
order to measure along + z (and — z). The corresponding 
operators are given by 

PP> 1= 0, a2 = 0) = |I)(I| 
Pl z (ai=n,a 2 = 0) = |II)(II|, (21) 

The results of the projective measurement are plotted 
versus different angles S of the spin analysis, which is 
depicted in Fig. [3] Complementary oscillations were ob- 
tained due to the spin flip in path |II). These curves are 
insensitive to the geometric phase 7, due to the lack of 
superposition with a referential sub-beam. 

The angle a' is set by a superposition of equal portions 
of 1 1) and I II), represented on the equator of the Bloch 
sphere. The interferograms are achieved by a rotation 
of the phase shifter plate, associated with a variation 
of the path phase x, repeated at different values of the 
spin analysis direction S. The projective measurement for 
a'i = 7r/2, a' 2 = corresponds to a phase shifter position 
of x— (and a' 1 ± = a[ + tt = 3tt/2, a' 2 — to x = "")• 



1 




S(rad) 

FIG. 3: Typical intensity modulations obtained by inserting a 
beam block, being the projections on the ± z direction of the 
path measurement, denoted as P£ z : (cn —0,a2 =0) and 
: (cti = 7r,a2 =0). The oscillations remain the same 
when altering the geometric phase. 
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FIG. 4: Typical interference patterns of the O-beam (a[ — tt/2) for S = 0, tt/8, tt/4, tt/2, 37r/4, being the direction of the spin 
analysis, and geometric phase 7 = (left) and 7 = n/6 (right side). Intensities at the path phase \ = an d X — tt are extracted 
from least square fits of the oscillations. The resulting curves (rear panel) represent the projections to the ± x direction of 
the path subspace, denoted as PS_ X : (a'i = tt/2, a'2 = 0) and PY_ X : (a'f 1 = 37r/2,a2 = 0). The shift of the oscillations (see for 
instance 5 — tt/ '2), due to the geometric phase 7, yields a lower contrast of the curves P^ x and PE. X - 



Projection operators read as 

PfM = = o) = i((|i) + |ii))«i|+(ii|)) 

(22) 

PlM 1 - = = 0) = - |n»((i|-(ii|)). ~ 

The interferogram obtained for 7 = and 8 = tt/2, 
in Fig. QJ is utilized to determine the zero point of the 
path phase x, which defines the + x - direction (a'i = 
7r/2, a' 2 = 0) for the path measurement. 

In order to obtain phase shifter scans of higher accu- 
racy, scans over two periods were recorded (see Fig.|4]) 
and the values for x = and tt are extracted from the 
data by least square fits. These extracted points, mark- 
ing the ± x-direction of the path measurement, are plot- 
ted versus different angles of 8, as shown in Fig.0] rear 
diagram. 

All phase shifter scans were repeated for different an- 
gles 8 for the spin analysis from 8=0 to 8 = tt in steps 
of 7r/8, and for several geometric phases 7 (steps of 7r/6, 
and beginning form 7 = 7r steps of 7r/4 ), depicted in 
the rear panel of Fig. Q] for five selected settings of 8 
(8 = 0, 7r/8, 7t/4, 7t/2, 37r/4) and two geometric phases 
( 7 = 0,7r/6). 



2. Azimuthal Angle Adjustment 

Here the Bell angles (polar angles) remain fixed at the 
usual values and are set at 8 for the projective spin mea- 
surement, and by the beam block (and fixed phase shifter 
positions) for the path measurement. The angle between 
the measurement planes is adjusted by one azimuthal an- 
gle (a' 2 ), which is deduced by phase shifter (x) scans. 

For the spin measurement the directions are fixed and 



given by given by (3: /3\ = 7r/4, /3 2 = and j3': j3[ = 
3tt/4, 0' 2 = (together with (3^ = -3tt/4, f3' 1 ± = -tt/4, 
see Fig.Q] for Bloch description and Fig. [5] for measured 
interference patterns). For the projective path measure- 
ment the fixed directions read as ol\ — (af" = tt, see 
Fig. [3] for measurements with beam block), and a[ = tt/2 
(a'i = 3tt/2). Phase shifter (x) scans are performed in 
order to determine a' 2 , which is depicted in Fig. [5] for two 



y- y=n/6 




Q J 1 1 1 1 1 1 1 1 1 — 

<* ft 21C -ft0il2ft 



Path Phase, x(rad) 

FIG. 5: Typical interference patterns of the O-beam (a{ = 
tt/2) for 5 = 7r/4 = /3i, 5 = 3tt/4 = f3[, 5 = -3tt/4 = Pi 
and 8 = — 7r/4 = /3( x (P2 = P2 = 0) and geometric phase 
7 = (left) and 7 = n/6 (right). Phase shifter scans \ are 
performed for a forthcoming determination of a' 2 - 



7 



values of the geometric phase: 7 = and 7 = tt/6. One 
can see a shift of the oscillations due to the geometric 
phase 7. 



F. Data Analysis and Experimental Results 

1. Polar Angle Adjustment 

Using least square fits from the polar angle adjustment 
measurement curves in Fig. 3] and Fig. [31 together with 
Eq.fJ^D]) the S- value is calculated as a function of the 
parameters (3\ and j3[ which is plotted in Fig. [6] for 7 = 0, 
7 = 7r/2 and 7 = 7r ( 7 = and 7 = 7r are chosen since the 
fringe displacement is maximal for these two settings and 
7 = 7r/2 illustrates the increase of S to a value of 2). The 
local maximum of the surface is determined numerically. 
The settings for f3\ and (3[, yielding a maximal S- value, 
are compared with the predicted values for 0i and f3[ 
from Eqs. p3ap and (|13b[) , respectively. 




FIG. 6: S- values for different settings of geometric phases 
7, derived from the least square fits of the projective mea- 
surements along ± z (beam block Fig. [3} and ± x (varying \ 
Fig. [4} direction for the path measurement (ai = 0, a'i = n/2 
and «2 = a 2 — 0) using Eq. (|ll[l . f3i and f3[ represent the di- 
rection of the spin analysis, which are changed systematically 
by a variation of S while 02 and (j'2 remain constant at the 
value zero (polar adjustment). The position of the maximum 
is determined numerically for different settings the geomet- 
ric phases 7 (here for example 7 = where Smax = 2v2, 
7 = tv/2 with Smax = 2, and 7 = tt, with Smax = 2y2). 
The /3i and /3[ values result as predicted in Eqs. (|13a[) and 
(USB for Smax. 
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FIG. 7: (a) Polar adjusted S-values versus geometric phase 
7 with adapted Bell angles (j3i and f3[) according to the ge- 
ometric phase 7. (b) the corresponding modified Bell angles 
are plotted versus the geometric phase 7. 



The resulting S values, derived by using the adjusted 
Bell angles fit and {3[, are plotted in Fig.[7|(a) versus 
the geometric phase 7. The theoretical predictions from 
Eq. (fTT|) depicted as solid (color: green) line are evidently 
reproduced. The maximal S decreases from 7=0 to 7 = 
7r/2 where the boundary of the CHSH inequality 5 = 2 
is reached, followed by an increase to the familiar value 

s = 2V2. 

In Fig.[7](b) the deduced f}\ and (3[ values are plotted 
versus the geometric phase 7. j3i and j3[ follow the theo- 
retical behavior (solid line - color: green line), predicted 
by Eqs. (|I3ap and (|13b[) . One can see a peak for /3i (and 
a dip for f3[) at 7 = tt. 



2. Azimuthal Angle Adjustment 

In Fig. [5] we depict selected S values calculated from 
least square fits of the azimuthal angle adjustment mea- 
surements, where (5\ — it/A, [3[ — 37r/4, (3^ = 57r/4, 
Pi 1 - = -7r/4 and a[ = tt/2 (see Fig. [5]) and ct\ = 0, 
aj~ = tt (Fig. [3]) versus geometric phase 7. A simple shift 
of the oscillation of the S-value is observed due to the 
geometric phase (see Fig. [5] front panel) . The maximum 
S-value of 2y2 is always found for a' 2 — 7, which is 
indicated in the rear panel of Fig. [8] The complete mea- 
surement set of ^-values versus the geometric phase 7 
is plotted in Fig.\§K&)-S-value azimuthal adjusted. If no 
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adjustment is applied to ct2, which means a' 2 is always 
kept constant at a 2 — 0, S approaches zero at 7 = 7r 
and returns to the maximum value 2y/2 at 7 = 2tt (see 
Fig.^l(a)-S-value without adjustment). 

Figure [5](b) shows adjusted a' 2 versus the geometric 
phase 7: It is clearly seen, that adjusted a' 2 fulfills the 
theory condition (solid line - color: green line) namely a 
linear dependency as expressed in Eq. (|15ap . 



IV. DISCUSSION 

If no corrections are applied to the Bell angles the S- 
value decreases from 2^/2 at 7 = to zero at 7 = -k 
and regains the value of 2v2 at 7 = 2tt (Fig.[9](a), S- 
value without adjustment). Keeping the azimuthal angles 
fixed, an appropriate adjustment of the polar Bell angles 
determined by the geometric phase = arctan(cos7)) , 
yields a sinusoidal oscillation of the S-value (2 < S < 
2\/2, with period n, see Fig.JTJa)). Finally, the maxi- 
mum S-value of 2 a/2 can be observed, for all values of 
the geometric phase 7, if the difference of the azimuthal 
angles (angle between the analysis planes) equals the ge- 
ometric phase (a' 2 = 7), while the polar Bell angles re- 
main unchanged at the typical values for testing of a BI 
(Fig.[9ja), S-value azimuthal adjusted). 

In our experiment, unlike the proposed setup in [48| . 




FIG. 8: S- values derived from least square fits of the projec- 
tive spin and path measurements for /3i = 7r/4, f3[ = 37r/4, 
Pt = -3tt/4, Pi 1 - = -tt/4 and a[ = tt/2, a'i X = 3tt/2 (see 
Fig-O, and a\ = 0, — ir (see Fig.[3l, versus geometric 
phase 7. The maximum S-value of 2\/2 is always found for 
a 2 = 7 as predicted in Eq. (|15a[l (azimuthal adjustment). If 
no corrections are applied to the Bell angles (a' 2 = 0) S ap- 
proaches zero at 7 = n. 




Geometric Phase y(rad) 

FIG. 9: (a) Azimuthal adjusted S'-values versus geometric 
phase 7 with balanced Bell angle (a'2) according to the geo- 
metric phase 7, and without corrections, (b) the correspond- 
ing modified Bell angle is plotted versus the geometric phase 
7- 



the geometric phase is not acquired solely in one arm 
of the interferometer but by two successive rf flippers. 
Hence, an alternative approach towards generation of 
geometric phase is introduced here: The effect of the 
rf-flipper inside the interferometer is described by the 
unitary operator U(4>i), which induces a spinor rotation 

from Iff) to |^>, denoted as U(<kM) = 14)- The rota- 
tion axis encloses an angle <j>i with the x-direction, and 
is determined by the oscillating magnetic field B^- 1 ' = 
B^f cos(ujt + <j)i) ■ y. Without loss of generality one can 
insert a unity operator, given by 11 = ffi (4>o)U (4>o) > yield- 
ing 

7 

UfaM) = U(<t>i) UHMU(£o) \1t) (23) 

1 

where U(4>o) can be interpreted as a rotation from \^\) 
to with the x-direction being the rotation axis, and 
U'(4>o) describes a rotation about the same axis back 
to the initial state |f|~). Consequently, U '((pi)W (<f>o) can 
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be identified to induce the geometric phase 7, along the 
reversed evolution path characterized by 4>o, followed by 
another path determined by <pj. 

Due to the inherent phase instability of the neutron in- 
terferometer, it is necessary to perform a reference mea- 
surement for each setting of 7 and 6. This is achieved by 
turning off the rf-flipper inside the interferometer, yield- 
ing a reference interferogram. The oscillations plotted in 
Fig. [4] and Fig. [5] are normalized, by the contrast of the 
reference measurement, and the phase of the reference in- 
terferogram is taken into account (relative phase between 
the oscillations). 

At this point it should be noted that the average con- 
trast of ^50% (obtained for 6 = 7r/2 with maximum 
intensity of ~ 25 neutrons/sec.) is below the threshold of 
70.7%, required to observe a violation of a BI. 

Violation of a BI, for a spin-path entanglement in neu- 
tron interferometry, has already been reported in , the 
argument here is the influence of the geometric phase on 
the »S-value. Consequently a normalization as performed 
does not influence the validity of the results presented 
here. 

Finally we want to discuss some systematic errors in 
our experiment, in particular in the state preparation and 
in the projective spin measurement. Under ideal condi- 
tions no interference fringes should be obtained in the 
H-beam, due to orthogonal spin states in the interfer- 
ing sub-beams. Nevertheless we have observed intensity 
modulations with a contrast of a few per cent. This in- 
dicates, that the state preparation (rf flipper) was not 
perfect in some sense. The expectation values for the 
joint measurements Eqs.([7])- (fTTj) can be deduced for an 
arbitrary (spin) state, in the path of the IFM where the 
rf-flipper is located, 

|*Me„.(7)) - 4(|I) ® \1t) +e iX |H) ® e l7 (sin^) 

V2 V 2 (24) 

+ cos^»). 

Here 9 is determined by the fringe contrast in the H- 



beam. These systematic deviations from the theoretical 
initial state have been taken into account in the calcula- 
tion of the final S value. 

The asymmetry in the curve of the projective mea- 
surement along the ± x direction of the path measure- 
ment, denoted as P+ x : (a[ — n/2,a' 2 = 0) and P^_ x : 
(a 1 ^- = 3tt/2, a' 2 = 0) in Fig.|4]is considered to result from 
a misalignment of the static magnetic fields, at the posi- 
tion of the coil, such as the stray field of the first guide 
field, the second guide field and the two fields (x, z)- 
direction produces by the coil itself. 



V. CONCLUSION 

We have demonstrated a technique to balance the in- 
fluence of the geometric phase, generated by one subspace 
of the system, considering a BI. This is achieved by an 
appropriate adjustment of the polar Bell angles (keeping 
the measurement planes fixed) or one azimuthal angle 
(keeping the polar Bell angles at the well-known values), 
determined by a laborious measurement procedure. It is 
demonstrated in particular, that a geometric phase in one 
subspace does not lead to a loss of entanglement, deter- 
mined by a violation of a BI. The experimental data are 
in goo d agrement with theoretical predictions presented 
in [48|, demonstrating the correctness of the procedure 
as a matter of principle. 
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